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1 Introduction 



Consider the linear equation 

At = {A,^}, (1.1) 

where {A,tp} = ippA^ — ipxAp. Here the unknown function A{x,t,p) depends on the "spectral" 
parameter p and a function ip = tlj{U{x,t),p). Following [Ij, we call (11. ip the Vlasov-type 
equation generated by ip{U,p). Equation (11.11) is also known as the dispersionless Lax equation. 

It is easy to check that a partial hodograph transformation A{x,t,p) —>■ p{x,t,A) reduces 
(11. II) to the following conservative form 

Pt = tlj{U,p),. (1.2) 

Here A plays a role of parameter. 

For some functions ip Vlasov-type equations are closely related to integrable hydrodynamic 
chains [21 13] . A hydrodynamic chain associated with Vlasov-type equation can be derived by 
expanding A{x,t,p) at a singular point of the function ip. In such a case, formula (11. 2p yields 
conservation laws for the hydrodynamic chain. 

Example 1 {The Benney chain) [U O El [1] . The Vlasov equation (or collisionless Boltzmann 
equation) has the form 

At + Apf/, - pA, = 0, (1.3) 

2 

where ip = ^^ + U. Substituting the expansion 

A" 

A = p + — + — + — + — + ... 1.4 

p p^ p-^ p^ 

into (II. 3p . one derives the famous Benney hydrodynamic chain 

A'l = A'+' + kA'-'Al k = 0,1,2,..., 

where A° = U. Let 

— A _ ^ — ^ — 1^1 — 1^1 — 

^~ X~A2"A3"A4"- 

be the inverse series for (II. 4p . Functions Hi can be easily calculated: Ho = A^, Hi = A^ , H2 = 
A^ + [A^f, . . .. The formula 

--(^-). 

generates infinitely many conservation laws for the Benney chain: 

fc-i 
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Example 2 {The Kupershmidt chain) [TJIH]. The Kupershmidt hydrodynamic chain 

= 5^+1 + + (3kB^Bl, = 0, 1, 2, ... 

is connected to the Vlasov-type equation fll.ll) . where 

by the expansion 

, . B^ B^ 

where B^ = U. Conservation laws for this chain can be calculated in the same way as in Example 
1. 

The hydrodynamic chains described in Examples 1, 2 admit infinitely many hydrodynamic 
reductions [9|,[T0]. The corresponding Vlasov-type equations admit the same reductions. 

The following "integrable" functions ip were found in [2]: 

Case 1. ijj{U,p) = U + W{p), where 

W" = CiW^ + C2W' + C3 

and 

Case 2. ijj(U,p) = UW{p), where 



^-(^CiW'^ + C2W' + Cs 



W" 

W 

Here Cj are arbitrary parameters. The Benney chain corresponds to Case 1 with W{p) = p'^/2. 

In this paper we describe all possible "integrable" functions ilj{U,p) using the method of 
hydrodynamic reductions. The existence of hydrodynamic reductions have been proposed as a 
definition of integrability for dispersionless multi-dimensional equations in [9]. We apply this 
approach for Vlasov-type equations. 



2 Hydrodynamic reductions 

Suppose there exists a semi-Hamiltonian [TT] hydrodynamic-type system 

ri = v\r)ri t = l,2,...,N, (2.5) 

and functions U = u{r) and A = \{r,p) such that these functions satisfy f ll.ll) for any solution 
r(x, t) of system (12.51) . Then (12.51) is called a hydrodynamic reduction for the Vlasov-type equa- 
tion (11.11) . The partial hodograph transformation X{r,p) p{r, A) leads to the corresponding 
hydrodynamic reduction of (11.21) . 
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Substituting A = \{r,p) and U = u{r) in (11. ip and (11.21) . we obtain the equations 

At = tppXx - tpx>^p (2.6) 

and 

Pt = i){u,p)x. (2.7) 
Calculating the derivatives by virtue of (12. 5p . we obtain from (12. 6p that 

^ {\pi)u diU + [v\r) - i)p\di)^ ri = 0, 

where we use the notation di = d/dr^. Since r is arbitrary solution of system (12. 5p . we get 

XpiJudiU = [i:p~v'{r)]d,X, z = l,2,...,N. (2.8) 

Let us determine functions p*(r), i = 1,2, ...,N as solutions of the equations 

v\r)=ijp\p^p,. (2.9) 

Then (12.80 implies that the equation Ap = has solutions (pairwise distinct in the generic 
case), i.e. 

Ap|p=p. =0, t = l,2,...,N. (2.10) 

Without loss of generality we can fix the Riemann invariants r* of the system (12. 5p in such a 
way that 

r' = X\p=pi. (2.11) 
Indeed, if we substitute p = into equation (12.61) . then (12.61) . (I2.10p imply 

(A|p=pi)t = {}Pp\p=pi){X\p^pi)x- 

This means (see (12. 9p ) that \\p=pi satisfies (12. 5p and therefore \\p=pi = Riir'^) for some functions 
Ri. According to (I2.10p . the branch points of the Riemann surface determined by the equation 
A = A(r,p), are nothing but the Riemann invariants of system (12.51) . This fact is well-known 
for hydrodynamic-type systems that produced by the Whitham averaging procedure applied to 
multi-phase solutions of both integrable continuous dispersive equations and integrable discrete 
equations (see references in [3]). 

Substituting functions p(r. A), u{y) in (12.71) . we obtain 

= ^ I - _ ^ • (2-12) 
If we fix A = r^, k ^ i, then (12.111) implies p = p^ and we obtain 



Let us introduce the following notation: 

f. _ f., _ '^u\p=pk . 

J I /I I 1 J IK ,1 ,1 ; I. -/- rv. 

Wp\p=pi Wp Wp\p=p^ Wp\p=p'' 

The compatibility conditions dk{dip) = di{dkp), i ^ k are equivalent to the equations 

^2 fik dpkfk - fki dpifi + duifk - fi) + fi dpfk - fk dpfi /o 1 ^\ 

OikU = -^—j diudku. (2.14) 

Ji Ik 

Equations (12.131) . (12.141) form a system of equations named in [2] the generalized Gibhons-Tsarev 
system (cf. [lOj). 

Since u does not depend on p, ( ]2.14p implies the following functional equation 

rj ( fik dpkfk — fki dp^fi + du{fk — fi) + fi dpfk — fk dpfi ^^ _ (o -[^\ 

dp^ ^3;^ ) - (2.15) 

for the function ilj{u,p). In the next sections we study this functional equation and found its 
general solution. The general solution is expressed in terms of a pair of arbitrary solutions of 
the standard hypergeometric equation 

u{u - 1) y{u)" +[{a + p + l)u--f] y{u)' + a[3y{u) = 0. (2.16) 

Note that the general solution of the Chazy equation, which appears in the classification paper 
[12j . can also be parameterized by a pair of hypergeometric functions. Our solution ip{u,p) is 
a generalization of the solution h{C,,u) found in |T3] (see Example 4). 



3 Particular solutions 

In Sections 4, 5 we solve the functional equation ( 12.15^ in terms of quadratures of hypergeo- 
metric functions. In this section we consider some particular cases, where the result can be 
written more explicitly. 

Computing the numerator of the left hand side of (I2.15P and expanding it at = p, = p, 
we obtain that the vanishing of the coefficients up to 8-th degree in the corresponding Taylor 
series is equivalent to the following system of PDEs for the function %p{u^p): 

^^Ippi^l - ^tPpp'iPppptPpppp'iPl + i^lpi^pppppi^l - '^i^ppi^lppi^li^pu + 2i)l^ii)ppppijlii)pu+ 

(3.17) 

6 i^lpi^pppi^li^ppu - ^^pp-iplippppu - Q^pp^pu-ipuu + Q^ppipu^puu = o, 



'^'^ppp'^v^Pppu '4^pp'^pppp'4^u'4^PPU 'i'^pp'^ppp'^u'^^PPPU ~^ '^pp'^u'^PPPPU '^''Ppp'^ppp'^Pu'^pu'^v 

6 i^ppi'pui'uu + ^i^ppi'ui'ppui'uu + ^i^lpippppi^li^puu + Qi^ppi'ui'pui'puu - ^i^lpipli^ppuu = o. 



(3.18) 
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and 

'^'^ppp'^u'^pu'^uu ~l~ '^pp'^^pppp'^u^pu'^uu ^'^pp'4^ppp'^u'4^pu'4^uu ^'^pp'^pu^uu'^ 

3 i^ppi^pppi^li^ppui^uu + G'iPlpi^ui^pui^ppui^uu - ip^piplippppui^uu + Sippppiplippuu- 

(3.19) 

'4'pp'4^pppp'4^u'^puu ~\~ ^ '^pp'^ppp'4^u^pu'4^puu ~\~ ^ '^pp'4^u'4^pu4^puu 'i'^pp'4^u^ppu'4^puu 

3 tppplpppplp^lpppuu ''Ppp'^Pu'^pu'^Pppuu ~l" '^pp'^u'Ppppu-u 0. 

This system and the functional equation fl2.15p are invariant under any transformation of the 
form u —>■ f{u). They are also invariant with respect to the following symmetry group: 

"ip C2ip + cip + Co, p ^ k2ip + hp + ko. (3.20) 

These symmetries are associated with linear transformations of independent variables in (11.11) . 
(11.21) . Notice that any function ip = W{u) +pV{u) satisfies equation (12.151) and system (I3.17p - 

(EUD. 

The integrable cases ilj{u,p) = u + W{p) and tp{u,p) = uW{p) described in Introduction 
can be found directly from (I3.17p . We denote these CclSGS clS Case 1 and Case 2, respectively. 
We mention two more particular integrable cases: 

Case 3. 'ip{u,p) = W{u — p), where 

W" = CiW'^ + C2W'^ + CsW; 

Case 4. ip{u,p) = pu + W{p), where 

w" = ^(^ciw^ + C2W' + c^y 

Here q are arbitrary constants. The Kuperhsmidt chain belongs to Case 4. It is easy to see 
that Case 3 is connected to Case 1 by the transformation ip ^ p. In all these four particular 
cases system (I3.17p - (l3.19p is equivalent to some ordinary differential equation of the fifth order 
for the function W. 

System (I3.17p - (I3.19I) admits the substitution 

i;p = F{p,^), (3.21) 

corresponding to the factorization with respect to the symmetry group u f{u). As the result, 
one gets an overdetermined system of three PDEs for the function F{x,y). We do not present 
this system here because of its complexity. 

Case 1 corresponds to the case 1: F{x,y) = V{x), where V = W. Case 3 corresponds 
to 3: F{x,y) = V{y). Case 2 corresponds to 2: F{x,y) = yV{x), where V = W'/W. Case 4 
corresponds to 4: F{x, y) = y/x + V"(x), where V = W — W/x. 
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Another particular solution is determined by the function 



The corresponding equation (13.211) leads to equation (12 .Tp of the form 



Pt = dj,^u (p2 + ti) -^2, 
where are arbitrary constants. In this example the function F has the form 

(3.23) 

An investigation of this ansatz leads to the following result. The case, where or TV is a 
linear function, was considered above. Namely, the cases W = \,V = \ and V = y coincide 
with Cases 3, 1 and 2, respectively. The case W = x transforms to Case 2 by the substitution 
ip ^ p. If V" 7^ and W" ^ 0, then the following three classes of solutions (I3.23P exist: 

Case 5. VV" = 2V'' + dV + Ca, WW"2W'^ + CiW + Ca, 

Case 6. VV" = V'' + ciV + C2, WW" = W'' + ciW' + ca. 

Case 7. VV" = -V'' + ciV + C2, WW" = -W'' + dW' + C2. 

For all these cases the generic solution depends on 6 arbitrary parameters. Function (I3.22p 
belongs to Case 5 with ci = —3, C2 = 1. 



4 General solution 

Let us expand the left hand side of (I2.15P into the Taylor series at = p. Denote by ^(p*, p) the 
first nontrivial coefficient of this expansion. For fixed p consider 5* = 0, 5*^ = as a system of 
linear algebraic equations with respect to derivatives ipui'^^up^- Ks determinant does not vanish 
if ipui^pp 7^ 0. Solving this system, we obtain 

= (4.24) 

where Q is a polynomial with respect to ipp of degree not greater than 4 with coefficients 
depending on u only. Taking into account (14.241) . it is easy to extract from the equation Spi = 
that 



^Ip Qi^pV ^^-^^^ 

where i? is a polynomial of degree not greater than three. The compatibility condition of 
equations (14.240 and (I4.25P has the form 

Q'^Qxxx — RQQxx + RQl; — {R^ + RxQ)Qx + RQu — RxxQ^ + i^RR^ — Ru)Q = 0, (4.26) 
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where x = ifjp. 

Assume that the polynomial 

Q{x) = a{x - bi){x - b2){x - 63) (x - 64) 

has distinct roots and rewrite fl4.25p as 

'^PPP ^ ^1 _^ _^ ^4 ,^ 27-) 

i^pp i^p-h "' i^p-bi 

One can verify that the system consisting of equations fl3.17p -( l3TT9l) and (14.261) is equivalent to 
the following : 

1. The functions ki{u) are arbitrary constants such that A;i + ... + A;4 = 3. 

2. The functions a{u), bi{u) satisfy the following system of ODEs: 

b[ = a{l-ki)\{{bi-b,), 1 = 1,. ..A- (4.28) 

The function a{u) can be chosen arbitrarily due to the admissible transformations u — > s{u). 
Consider the double ratio 

_ (&l-&2)(&3-fe4) 
P ^b,-bs)ib2-b,)- 

Differentiating p by virtue of fl4.28p . it is easy to check that p' 7^ 0. Let us change u in such a 
way that p = u. This means that we choose 

1 1 



(&2-&3)(&l-&4) {b^-b2){b^-b,y 

One can verify that the formulae 

Z2 + uy2 I. y2 . Z2 + y2 , ^2 

bi = ■ , 62 = — , 03 = ■ , 64 = — , 

Zl + uyi yi Zl + yi zi 

where (yi, Zi), (1/25 ^2) are two arbitrary solutions of the linear system 

, ki + k2 + k3-2 ki + k2 + k3-2 , I - ki 1 - fci - A;2 + ^2^ , , 

y= 1 y + 7 7^ z= -y + z, 4.29 

u — 1 u[u — 1) u — 1 uyu — 1) 

define a general solution of (14. 280 . Notice that if ki + k2 + k^ ^ 2, then 

u{u — 1) I 

z = -uy + - — y 

ki + k2 + k3-2 

and system (I4.29P is equivalent to the hypergeometric equation (I2.16p . where /ci = 1 + a — 7, 
^2 = 1 - a, A;3 = 7 - /?. 
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System (14.241) . fl4.25p can be reduced to quadratures by the following way. Let us determine 
a function 4>{u,p) as the solution of the over-determined system: 

-^)y'i . _ 0^1(0 ^^^^^^ 



(3{yi(j) + zi)' yi(p + zi 

It is easy to check that this system is consistent. Then the solution of the following system in 
involution 

(3{yi(j) + zi) yi4> + zi 

is a general solution of (12.151) . This fact can be verified by a direct calculation. It turns out 
that the expression under differentiating in the left hand side of (12.151) is equal to 

{I - u - P)a'^ + (1 - u - a)/?^ + u (4a/3 - a - /?) 

l)(a-/5)2 ' 

where a = (j){u,p^),P = (f){u,p^). 

Remark. The standard Wronskian formula for second order linear ODE implies that the 
expression y2y'i — yiy'2 from (I4.3ip equals Cu^^^^^~'^{u — i)'=2+fc3-2 f^j, some constant C. 



5 Degenerations 

In Section 4 we have considered the general case. This means that the polynomial Q has 
degree 4 and all its roots hi are distinct for the generic value u. In this section we consider 
degenerations. It is easy to see that the degree of the polynomial Q can be fixed by 4 with the 
help of transformations (13.201) . It turns out that in all cases the result can be parameterized 
by a pair of solutions of some degenerations of the hyper geometric equation. 

Degeneration 1. Suppose Q = a{x — bi)'^{x — b2){x — 63); then 



Q x — bi (x — 6i)2 X — 62 x — b^' 
where ki, ^2, ^3 are constants such that ki + k2 + = 3, and 

b[ = -a{b, - b2){b, - 63)/!, b', = a{b2 - b,f{b2 - b,){l - k^), 

h'^ = a{b^-b^f{b^-b2){l-k^), /( =a/i((6i-62)(6i-63)(2-A:i) + (62 + &3-26i)/i). (5.32) 
A general solution of system (15.321) can be parameterized in the following way: 

, z/i , ^1 u y^ + 

02 = — , 03 = — , oi = ■ , 

2/2 Z2 y2 + UZ2 
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where {yi,zi), (^2,-22) are two arbitrary solutions of the hnear system 

The functions a, fi are determined by system (15.321) . Note that system fl5.33p is equivalent to 
the Bessel equation 

//^ /I (^2-^1X^1+^2-2) • 



Degeneration 2. Suppose Q = a{x — bi)^{x — 62)^; then 



Q X — bi (x — 61)^ a; — 62 (x — 62)^' 
where ki, k2 are constants such that ki + k2 = 3, and 

b[ = -a{b, - b^fh, b', = -a(6i - b^Yh, 

f[ = a(62 - &i)/i ((&2 - &i)(A:2 - 1) + 2/i) , = a(6i - 62)/2 ((^i - 62)(fci - 1) + 2/2) . 

The general solution is given by 

bi = — , 02 = — , 

where {yi,zi), {y2,Z2) are two arbitrary solutions of the linear system 

' 3A:i-5 , 1 , 3(1 -fci) 

y=— y + uz, z=--y + — z. 5.34 

2u S 2u 

Notice that the function y{u) satisfies the following second order equation: 

I P 3(3A:,-5)(3A:,-7) 
\3 



Degeneration 3. Suppose Q = a{x — bi)^(x — 62); then 

R ki fi f2 ^ k2 



Q X — 61 (x — 61)^ (x — 61)'^ X — 62' 
where /ci, /c2 are constants, such that ki + k2 = 3, and 

fe'i = a{b2 - 6i)/2, 62 = a{b2 - 6i)=^(l - k2), 

f[ = a((fe2 - &i)(/i' - ^2/2) - 2/1/2), /2 = 2a/2((62 - h)f, - /s 
The general solution can be written in the following form 

, yi , y'l + uyi 

bi = — , 62 = , 

2/2 2/2 + uy2 
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where yi, 1/2 are arbitrary solutions of the Weber equation y" = + (1 — 2A;2) jy. The functions 
^5 fl, /2 are completely determined by the above system. 
Degeneration 4. Suppose Q = a{x — 6)^; then 



R h . h . f2 ^ h 



Q X — b (x — b)^ (x — b)^ {x — b)^ ' 
where ki = 3, and 

b' = -afs, il = -3ahh, r, = -2aUi + hh), f[ = a{h-2hh). 
Eliminating /a, /2, /i from this system and choosing a = —1/6'^, one obtains the equation 

b'" _ 3 6"2 _ _ 

Its general solution can be written in the form b = yi/y2, where yi,y2 are arbitrary solutions 
of the Airy equation y" = uy. 

Deeper degenerations can be obtained by the restriction that the polynomial Q possesses 
one or several constant roots. In this case it is convenient to make one of these roots the 
infinity using transformation (13.20p . Consider, for instance, Degeneration 4 under assumption 
b = const. Choosing the normalization a = 1, one obtains /s = and 

= -2fl f[ = -2/1/2. 

The simplest solution /2 = /i = of this system corresponds (for b = 00) to the Benney 
pseudo-potential from Example 1. The solution /2 = 0, /i 7^ implies ip{u,p) = up + plog (p), 
which coincides with Case 4 from Section 3 for ci = C2 = 0, C3 = 1. Finally, if /2 7^ 0, one can 
obtain (up to the equivalence) 



In this case the solution cannot be expressed in terms of elementary functions. 

For general system fl4.28p the fact that some roots bi are constant, is equivalent to the 
equality fcj = 1 for corresponding values k^. Solutions of system 04.281) for such degenerations 
can be extracted from fl4.30l) . fl4.3ip . We omit the explicit formulae for such cases and the 
analysis of the case of constant bi for Degenerations 1-3. 

Let us describe particular solutions from Section 3 in the context of Sections 4, 5. It turns 
out that Case 5 with cf 7^ 8c2 is equivalent to fl4.28p . where bi = —62, &3 = — &4, ^2 = ki 
ki = k^. Namely, 

3C3 - Ci 3C3 + Ci 

ki = — , fca = — . 

4C3 4C3 

for C2 = cj/S - 4/8. 
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Case 6 with 7^ 4c2 is equivalent to (14.281) . where two roots bi are constant. The right 
hand side of (14.251) reduces to the form 

1 h k2 
- + ^ + 



(5.35) 



X X — hi X — 1)2 

where x = ipp, 

J C3- Ci C3 + C1 

K^ = , fco = . 

2C3 ' 2C3 

for C2 = Ci/4 - C3/4. 

In Case 7 with 7^ — 4c2 we have Q{x) = ax{x — &i)(x — 62) and 

R ki k2 
Q X — hi X — 62 ' 

where 

3C3 - Ci 3C3 + Ci 

ki = , ko = 

2C3 ' 2C3 

for C2 = — c^/4 + C3/4. Notice that in this case fci + ^2 = 3. It is easy to verify that if the right 
hand side of ( 14.25^ has the form (]5.35p . then the constants ki,k2 can be arbitrary, and 61,62 
satisfy the system 

b[=a{hi-h2){l-ki), h'^ = a{h2-hi){l-k2). (5.36) 

In the case, when ki + k2 = 2, system (I5.36P possesses the solution a = l;6j = u + ti, 
where ki = 1 + l/(ti — ^2), ^2 = 1 + 1/(^2 — ^i)- It corresponds to Case 4 from Section 3 with 
Ci 7^ 0, C2 7^ 4C1C3. For the Kupershmidt chain (see Example 2), 

R ki 



Q X — hi 

a = 1 and hi = u + ti. 

For Case 2 with ci 7^ 0, C2 7^ 4ciC3 we have (15.351) . (I5.36P and ki + k2 = 1. Under the 
latter condition, (I5.36P possesses the solution a = l/u; hi = tiU, where ki = ^2/(^2 — ^1)5^2 = 
ti/iti-t2). 

Case 1 corresponds to a constant solution of system (I5.36p . which exists for ki = k2 = 1. It 
was already mentioned in Section 3 that Case 3 is equivalent to Case 1. 



6 Conclusion 



We apply the method of hydrodynamic reductions to classify integrable Vlasov-type equations 
of the form (11.11) . (11.21) . In this paper the simplest case of one function U{x,t) is completely 
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analyzed. In the next paper we are going to solve a more complicated problem of classifica- 
tion of integrable Vlasov-type equations in the case of two functions Ui{x,t), U2{x,t). It turns 
out that there exist several essentially different classes of integrable functions iIj{Ui,U2,p)- 
One of such classes corresponds to two-component (2+l)-dimensional hydrodynamic-type sys- 
tems. This class was constructed in the paper [H]. Note that examples of integrable functions 
il^{Ui, . . . , Un,p) appeared earlier in other papers, where other approaches were used. In partic- 
ular, functions ip, associated with algebraic curves of an arbitrary genus, were constructed in 
[T3] . An integrable function ip{Ui, . . . , Un,p) was constructed from any n-dimensional Frobenius 
manifold in |16j . 

As it was shown in [10], equation f l2.12p in the case of the Benney chain (see Example 1) 
is nothing but the Loewner equation well known in the theory of conformal mappings. The 
results obtained in our paper can be of interest in connection with the so-called Laplacian 
growth problem (see [TTj and references therein). Moreover, each integrable case leads to an 
integrable hydrodynamic chain similar to the Benney chain (see Example 1.) If the range of 
the discrete variable k is the set of all integers, the corresponding hydrodynamic chains can be 
constructed rather easily. However, the problem of a "right" truncation of such chains to the 
set of non-negative values of k is not trivial. We are going to write a separate paper on the 
subject. 
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